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ABSTRACT  ■  .* 

Measurements  of  three  wave  surface  profiles  in  a  line 
perpendicular  to  the  beach  extending  across  the  breaker  line 
and  a  set  of  horizontal  and  vertical  water  particle  velocities 
at  the  breaker  line  were  taken  in  a  mild  surf  and  spilling 
type  breakers.   The  application  of  a  linear  wave  theory 
transfer  function  to  the  surface  profile  spectrum  resulted 
in  a  calculation  of  the  horizontal  water  particle  velocity 
spectrum  which  was  13  percent  lower  than  the  measured 
spectrum  peak  and  showed  a  phase  five  degrees  ahead  of  the 
measured  horizontal  water  particle  velocity.   Celerity 
calculations  between  adjacent  wave  profiles  over  the  coherent 
frequency  range  of  the  spectra  showed  that  the  wave  phase 
speed  was  not  frequency  dependent  and  was  predictable  by  the 
linear  theory  shallow  water  approximation.   The  finding  that 
a  dispersive  linear  wave  system  predicts  the  wave  kinematics, 
while  a  non-dispersive  system  was  observed,  suggested  a 
contradiction  and  was  discussed.   Capacitance  wave  gages 
and  electromagnetic  flow  meters  were  found  to  be  excellent 
instruments  for  use  in  the  surf  zone. 
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I.   INTRODUCTION 

A.  OBJECTIVES  OF  THE  EXPERIMENT 

The  kinematics  of  wave  forms  and  water  particle  vel- 
ocities in  the  surf  and  breaker  zones  has  been  an  enigma 
to  science  for  many  years.   Laboratory  studies  have  estab- 
lished the  fact  that  the  linear  solution  to  the  wave 
equations,  which  applies  well  to  waves  of  moderate  height 
in  deep  water,  cannot  fully  describe  monochromatic  waves 
and  associated  water  particle  velocities  during  the  break- 
ing of  waves.   This  experiment  was  designed  to  measure 
surface  profiles  and  water  velocities  in  the  breaker  zone 
and  compare  the  results  with  present  theory.   As  a  secondary 
objective  the  accuracy  and  usefulness  of  capacitance  wave 
gages  and  electromagnetic  flow  meters  in  the  surf  environ- 
ment were  to  be  evaluated.   The  experiment  also  added  to 
a  data  bank  of  measurements  taken  in  the  surf  zone  for 
studies  being  conducted  at  the  Naval  Postgraduate  School. 

It  is  expected  that  the  data  obtained  will  ultimately 
contribute  to  the  establishment  of  a  theory  which  will 
accurately  predict  various  wave  parameters  from  one  set  of 
measurements  taken  in  the  breaker  zone. 

B.  APPROACH  TO  THE  PROBLEM 

Measurements  of  water  surface  profiles  at  four  locations 
in  a  line  perpendicular  to  the  beach  and  two  sets  of 
horizontal  and  vertical  water  particle  velocities  at  the 
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breaker  point  were  attempted  during  this  research.   An 
energy-density  spectrum  was  ca'i'cui:ated  for  each  data  set 
and  coherence  and  phase  were  calculated  between  data  sets. 
A  wave  profile  spectrum  was  converted  to  horizonal  and 
vertical  velocity  spectra  using  linear  theory  transfer 
functions.   These  spectra  compared  with  the  measured  ve- 
locity spectra  taken  at  the  same  point.   Solitary  wave 
theory  was  also  compared  to  the  measured  results.   The 
assumption  that  the  processes  in  the  breaker  zone  are  sta- 
tionary and  ergotic  was  investigated. 

C.   PREVIOUS  WORK 

This  experiment  is  a  continuation  of  work  conducted  at 
the  Naval  Postgraduate  School  by  Steer  (1972^  and 
Richardson  (1973)  under  the  direction  of  Thornton  (1968, 
1972).   The  problem  of  taking  direct  measurements  in  the 
surf  zone  was  first  attacked  by  Inman  (1956)  who  infered 
water  particle  velocities  from  measured  drag  forces  under 
the  waves.   Miller  and  Zeigler  (1964)  measured  in  situ 
particle  motions  using  acoustic  and  electromagnetic  current 
meters.   Walker  (1969)  made  the  same  type  of  measurements 
using  the  propeller  type  flow  meters.   Wood  (1973)  measured 
waves  and  currents  in  the  surf  zone  using  movies  of  dye 
movement  and  capacitance  wave  gages.   Work  by  Dean  (1974) 
gives  an  analysis  of  various  wave  theories  applicable  in  the 
surf  zone. 
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II.   DATA  COLLECTION  AND  PREPARATION 
A.   EQUIPMENT 

Measurements  were  taken  using  two  electromagnetic  flow 
meters,  three  capacitance  wave  height  gages,  and  one 
resistance  wave  height  gage.   All  equipment  was  calibrated, 
in  the  laboratory,  prior  to  the  experiment. 

The  two  flow  meters  were  Marsh-McBirney  Model  721 
Electromagnetic  Current  Meters.   The  flow  meter  operation 
is  based  on  Faraday's  principal  of  electromagnetic  induc- 
tion.  Each  probe  measures  water  velocity  in  two  orthagonal 
directions,  horizontally  (shoreward)  and  vertically. 
Measurement  accuracy  was  determined  to  be  +0.02  m/s  during 
calibration.   A  carpenters  level  was  used  to  establish  axis 
alignment  in  the  field  with  an  estimated  error  of  +5 
degrees. 

The  flow  meters  were  calibrated  in  the  laboratory  using 
an  oscillating  platform  attached  to  an  eccentric  arm  50.72  cm 
long  driven  by  a  variable  speed  motor  using  the  method  of 
Steer  (1972)  and  Richardson  (1973).   The  speed  of  the  motor 
was  determined  using  a  stopwatch  to  measure  the  time  re- 
quired for  five  revolutions  and  then  making  a  comparison 
with  the  average  of  ten  maxima  voltage  outputs  measured  on  a 
strip  chart  recorder.   Because  the  platform  became  unstable 
at  velocities  above  1.3  m/s  (0.4  Hz),  calibration  stopped 
at  this  point.   The  flow  meter  calibration  plots  are  shown 
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in  Figure  (1) .   A  linear  regression  was  performed  on  the 
actual  velocity  versus  voltage  measurements  to  obtain  the 
calibration  factors  listed  in  Table  (1) . 

The  three  capacitance  wave  staffs  were  constructed  by 
the  author  using  the  design  of  McGoldrick  (1969) .   The 
2-1/2  m  gages  consisted  of  a  1.0  cm  diameter  stripped 
coaxial  cable  mounted  on  a  painted  wooden  frame  as  shown  in 
Figure  (2) .   The  capacitance  wave  gage  operates  on  the 
principal  that  a  change  of  capacitor  plate  dimensions 
changes  the  circuitry  voltage  output.   The  insulated  coaxial 
cable  and  the  sea  water  act  as  the  plates  and  the  cable 
insulation  as  the  dielectric.   As  the  wave  elevation  changes, 
the  capacitance  of  the  circuit  changes  and  the  voltage  output 
responds  linearly. 

The  frame  was  rather  heavily  constructed  and  attached 
close  to  the  tower  to  maintain  a  fairly  rigid  probe  wire  and 
to  reduce  water  resistance  on  unattached  parts.   When  assembled, 
the  frame  became  somewhat  unwieldy  in  the  surf  and  created 
some  resistance  to  flow  around  the  support  tower.   The  accuracy 
of  the  gages,  however,  was  not  believed  effected.   The 
electronic  gear,  including  power  source,  was  mounted  in  a 
waterproof  brass  container  on  top  of  the  frame  and  two  conductor 
cables  carried  the  output  signal  shoreward. 

The  wave  gages  were  calibrated  in  the  laboratory  to  a 
depth  of  1.7  m.   Accuracy  was  estimated  to  be  +0.005  m.   The 
calibration  plots  are  shown  in  Figure  (1)  and  the  results 
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TABLE  (1)  -  INSTRUMENT  DESIGNATIONS 
AND  CALIBRATION  INFORMATION 


TRACK  # 
TITLE 

INSTRUMENT 

1 

Outer 
Profile 

WS2 

2 

Center 
Profile 

WSl 

1 

3  . 

Inner 
Profile 

BAY 

4/5 

Horiz/Vert 
Velocity 

FM1X/FM1Z 

MANUFACTURER 
MODEL/SERIAL# 

Locally 
Made 

Locally 
Made 

Baylor 
13528R 

Marsh- 
McBirney 
721/95 

RANGES 

VOLTAGE  OUT 

MEASUREMENT 

3.40-8.29v 
0.0-2. 5m 

2.40-5.35V 
0. 0-2. 5m 

0.33-0.16v 
0.0-1. 52m 

+1.0v 

,'  +3 .  05m/sec 

HEIGHT  ABOVE 
BOTTOM 

0.74m 

0.33m 

0.24m 

0.53m 

LAB  CALCULATIONS1 
MULT  FACTOR 
ADD  FACTOR 
STD .  ERROR 

+0.511 
-1.753 
±0.  004m 

+0.846 
-2.029 
+0.007m 

-9.405 
+3.078 
+0.020m 

+10.908 
-0.223 
+0. 014m/sec 

PROCESSING  GAINS2 
A 
B 
C 

xO.177 

x50 
xO.1131 

xO.236 

x50 
xO.0848 

X3.536 

x50 
xO.0057 

xl.414 
xlOO 
xO.0071 

FINAL  CALIBRATION3 
MULT  FACTOR 
ADD  FACTOR 

i 

+0.058 
-0.995 

+0.072 
-1.699 

-0.053 
+3.313 

+0.024 
-0.068 

From  linear  regression  of  laboratory  calibrations  as 
shown  in  Figure  (1) . 

2  A  -  Instrument  to  analog  tape,  B  -  Analog  tape,  through 
filtering,  to  digital  tape,  and  C  -  1/ (AxB) ,  applied  to 
digitally  recorded  voltage  to  compute  instrument  output. 

3  Factors  applied  to  data  in  digital  form  (including 
instrument  height)  during  computer  analysis  to  obtain  actual 
surface  height  and  velocity  values. 
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of  a  depth  versus  voltage  output  linear  regression  are 
shown  in  the  calibration  factors  of  Table  (1) . 

The  resistance  wave  gage  was  a  Baylor  Wave  Measuring 
System  Model  #13528R  with  a  1.52  m  wave  height  capacity. 
The  gage  consisted  of  two  1.27  cm  stainless  steel  wire 
ropes  mounted  under  tension  on  a  tower  and  an  electronics 
box  located  ashore.   As  the  water  level  rose  on  the  cables, 
they  were  shorted  and  a  voltage  was  produced  which  was 
linearly  proportional  to  the  amount  of  exposed  cable. 

Calibration  consisted  of  simply  shorting  the  cables  at 
various  intervals  and  recording  the  output  voltages.   Cali- 
bration accuracy  was  +0.02  m  when  measurements  were  taken 
away  from  either  end  of  the  gage.   The  calibration  plot  is 
shown  in  Figure  (1)  and  the  calibration  factors  determined 
by  linear  regression  are  shown  in  Table  (1) . 

All  data  was  recorded  on  a  14  channel  Sangamo  Analog 
Recorder.   Voltage  inputs  were  reduced  to  a  range  of  -1.0 
to  +1.0  v  and  recorded  continuously.   A  strip  chart  recorder 
was  used  to  monitor  both  input  and  recorder  output  voltages 
during  the  experiment.   One  wave  profile  data  set  (WS3)  was 
lost  during  recording. 

B.   BEACH  SET  UP 

The  experiment  was  conducted  at  Del  Monte  Beach  at  the 
Naval  Postgraduate  School  Beach  Laboratory  on  Thursday, 
23  May  1974.   Data  was  recorded  from  12:50  to  14:30  PM  (PDT) 
(100  minutes).   High  tide  was  at  13:40  PM.   The  weather  was 
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clear,  with  some  gusty  onshore  wind  (less  than  10  knots)  and 
no  visible  sea.   The  swell  surf  was  mild  with  observed 
significant  wave  heights  of  between  0.5  and  0.75  m  and  periods 
of  approximately  8  to  10  seconds.   Breakers  were  predominantly 
spilling  with  an  occasional  plunging  type.   Refraction  along 
the  shore  consistantly  caused  waves  to  break  paralled  to  the 
beach. 

The  four  towers  used  were  steel  pipes  6.3  cm  in  diameter 
and  3.6m  high  with  a  1.0  m  diameter  base.   They  were  held  up 
with  0.9  cm  polypropylene  rope  tied  to  homemade  screw-type 
anchors.   A  typical  tower  arrangement  is  shown  in  Figure  (3). 
The  towers  were  placed  on  a  line  perpendicular  to  the  shore 
at  specified  distances.   The  towers  remained  visibly  stable 
during  measurements  and  it  was  felt  that  tower  movement  was 
not  appreciable.   The  polypropylene  lines  had  a  tendancy  to 
stretch  and  loosen  and  knots  did  not  hold  well.   The  use  of 
steel  guy  wires  would  have  reduced  line  slackness  problems. 
The  anchors  were  difficult  to  both  install  and  remove  and 
did  not  hold  unless  deeper  than  about  0.5  m.   The  experiment 
was  concluded  when  Tower  Two  became  unstable  and  fell  as  a 
knot  failed.   Flow  meter  number  two  (FM2)  was  mounted  too 
high  on  the  tower  and  was  occasionally  exposed,  causing 
false  readings.   Thus  this  data  could  not  be  used  for 
continuous  analysis.   Figure  (4)  is  a  diagram  of  the  beach 
profile  and  the  instrument  arrangement. 
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C.   DATA  PREPARATION 

Data  was  recorded  at  the  beach  on  eight  channels  at  a 
speed  of  1  7/8  in/sec.   The  data  was  played  back  at 
7  1/2  in/sec  and  sampled  at  100  Hz  (1  sample  per  0.01 
seconds)  or  25  Hz  real  time  (1  sample  per  0.04  seconds). 
The  analog  to  digital  conversion  was  performed  on  a  CDC/9300 
computer.   A  summary  of  tape  recorder  and  digitization  gains 
is  included  in  Table  (1) .   Playback  data  was  amplified 
once,  filtered,  amplified  again,  digitized,  and  recorded. 
A  matched  set  of  Krohn-Heit  filters  with  a  cutoff  frequency 
of  4  0  Hz  (10  Hz  real  time)  was  used  to  low  pass  filter  the 
analog  signals.   Filtering  was  conducted  during  this  phase 
to  reduce  the  possibility  of  60Hz  aliasing.   The  real  time 
sampling  interval  was  0.04  seconds  so  the  recorded  Niquist 
frequency  was  12.5  Hz.   Digitized  data  was  recorded  on  a 
nine  track  magnetic  tape  (labeled  BUB191)  in  one  file  of 
597  records,  each  record  consisting  of  2048  points  or  256 
data  points  for  each  of  eight  channels.   Figure  (5)  gives 
examples  of  the  recorded  data  taken  during  the  experiment. 
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III.   SPECTRAL  ANALYSIS 

A.   APPROACH  TO  DATA  ANALYSIS 

Continuous  time  series  records  of  all  measurements  were 
taken  for  100  minutes.   A  thirty  minute  sample  record  of 
each  data  set  was  spectrally  analyzed  on  an  IBM/360 
computer.   The  power  spectra  generated  described  the  com- 
position of  the  data  over  a  frequency  range  in  terms  of  its 
mean  square  value.   An  examination  of  a  power  spectrum 
indicates  the  regions  of  greatest  mean  wave  or  velocity 
amplitude  and  their  corresponding  frequencies  and  will  give 
a  general  relationship  of  energy  concentrations  across  the 
frequency  range.   A  comparison  of  two  spectra  ?jy  means  of 
a  cross-spectrum  and  subsequently  a  coherence  and  phase 
plot  with  frequency  will  indicate  the  regions  and  degree  of 
linear  relationship  and  phase  between  two  data  sets.   Co- 
herence can  be  thought  of  as  the  ratio   between  two  measure- 
ments of  the  square  of  the  gain  factors  for  a  constant 
parameter  linear  system.   A  coherence  near  1.0  indicates 
that  both  spectra  operate  under  the  same  gain  factor  or 
transform  function,  that  is,  their  relationship  is  highly 
linear.   A  coherence  near  0.0  indicates  the  two  examined 
spectra  are  not  related  or  are  related  only  by  a  non-linear 
gain  factor.   Phase  angle  measurements  indicate  the  time 
delay  between  two  data  set  spectra. 
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The  use  of  spectral  analysis  in  the  examination  of  data 
required  the  following  assumptions  be  made: 

1.  The  data  was  random.   In  other  words,  the  data  used 
represented  one  of  many  possible  results  that  could  have 
occured  or  was  one  physical  realization  of  the  random  data 
process. 

2.  The  random  process  was  stationary.   That  is,  the 
data  set  could  be  described  in  terms  of  certain  average 
values  which  are  true  for  all  data  sets  of  the  same  process. 

3.  The  stationary  random  process  was  ergotic  or  all 
properties  of  the  process  for  all  time  could  be  described 
by  performing  analysis  on  one  finite  data  set. 

Many  processes  found  in  nature  can  be  treated  as  being 
ergotic,  stationary,  and  random.   The  validity  of  this 
assumption  as  it  applied  to  this  experiment  will  be  examined 
in  Section  (IV, C)  of  this  paper. 

B.   COMPUTER  ANALYSIS  METHOD 

A  complete  flow  chart  of  data  acquisition  and  analysis 
is  presented  in  Figure  (6)  and  is  discussed  briefly  below. 
The  equations  used  are  given  in  Bendat  and  Piersol  (1971) . 
The  parameters  used  for  the  spectral  analysis  process  were 
chosen  to  balance  the  longest  record  which  the  computer 
could  reasonably  analyze  with  the  best  resolution  over  the 
frequency  range  of  interest.  The  frequency  band  of  0.0  to 
1.5  Hz  was  chosen  to  cover  the  region  of  gravity  waves 
(0.033  to  1.0  Hz).  Low  pass  filtering  was  used  twice  to 
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Figure     (6) 


minimize  60  Hz  aliasing.   Folding  of  60  Hz  noise  around  the 
Niquist  frequency  would  have  appeared  as  peaks  in  the 
spectra  at  0.625  and  0.9375  Hz.   No  such  peaks  showed  up 
during  analysis,  demonstrating  the  success  of  the  filtration 
employed. 

The  same  30  minutes  of  raw  data  was  analyzed  for  each 
of  the  five  channels.   During  each  program,  two  channels 
were  chosen,  digitally  filtered,  and  decimated  by  eight.  The 
data  was  low  pass  filtered  using  a  25  weight  filter  with  a 
cutoff  frequency  of  1.2  Hz  and  terminal  frequency  of  1.8  Hz. 
This  filter,  illustrated  in  Figure  (7),  is  given  in  Davidson 
(1970).   The  final  analyzed  time  step  was  thus  0.32  seconds 
with  a  Niquist  frequency  of  1.5625  Hz. 

The  data  points  were  then  converted  to  wave  height 
and/or  velocity  values  using  the  calibration  factors  of 
Table  (1) .   A  mean  value  was  calculated  for  both  data  sets 
and  the  data  was  linearly  detrended.   The  variance,  standard 
deviation,  and  average  period  were,  calculated.   The  average 
period  was  determined  by  finding  the  time  between  zero  axis 
upcrossings;  the  calculated  average  periods  appear  low 
using  this  method.   The  mean  values  for  each  data  set  are 
listed  in  Table  (2) . 

For  pairs  of  measured  parameters,  the  auto-covariance 
functions  were  calculated  using  a  five  percent  lag.   Each 
auto-covariance  was  smoothed  with  a  Parzen  window.   A 
Fourier  transform  was  then  applied  to  the  smoothed  auto- 
covariance  functions  and  the  two  power  density  spectra 
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TABLE  (2)  -  STATISTICAL  PARAMETERS 
CALCULATED  DURING  ANALYSIS 


TRACK 
TITLE 

1 

OUTER 
PROFILE 

2 

CENTER 
PROFILE 

3 

INNER 
PROFILE 

4 

HORIZONTAL 
VELOCITY 

5 

VERTICAL 
VELOCITY 

DESIGNATION 

WS2 

WS1 

BAY 

FM1X 

FM1X 

MEAN  DEPTH 
OR  VELOCITY 

(m  or  m/sec) 

1.982 

1.468 

1.428 

+0.018 

+0.048 

VARIANCE 
(m2  or  m2/sec2) 

0.016 

0.029 

0.032 

0.170 

0.006 

STANDARD 
DEVIATION 

(m  or  m/sec) 

0.127 

0.170 

0.178 

0.412 

0.077 

AVERAGE  PERIOD 
(seconds) 

5.53 

5.29 



4.20 

< 

6.20 

_-l 

3.10 
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determined.   Bandwidth  resolution  for  the  spectra  was 
0.00544  Hz. 

A  cross-covariance  function  between  data  sets  was  then 
computed,  smoothed  with  a  Parzen  window,  a  Fourier  trans- 
form applied,  and  the  cross-spectrum  computed.   The 
coherence  and  phase  relations  were  calculated  and  the  two 
spectra,  the  coherence,  and  the  phase  plotted. 
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IV.   RESULTS 

A.  INTRODUCTION 

The  spectra  generated  were  examined,  the  observed  low 
frequency  peaks  were  related  to  phenomenon  which  one  would 
expect  to  find  in  the  surf  zone.   The  stationarity  of  the 
data  sets  was  examined  in  order  to  determine  the  scope  of 
applicability  of  the  data. 

An  examination  of  theories  which  may  be  used  in  the 
surf  zone  was  conducted  using  two  approaches.   First, 
measured  velocity  spectra  were  compared  with  velocity  spectra 
calculated  from  the  center  surface  profile  using  a  linear 
wave  theory  transfer  function.   Solitary  wave  theory  was 
then  to  be  used  to  develop   a  discrete  line  spectrum  for  a 
similar  type  of  comparison.   However,  the  results  showed 
that: (1)  the  linear  theory  transfer  function  calculated  hori- 
zontal velocity  very  well  and  (2)  there  were  no  identifiable 
harmonics  in  the  measured  spectra  which  could  be  related 
to  such  a  discrete  spectrum;  this  indicated  that  solitary 
wave  theory  was  not  applicable  under  the  experiment 
conditions . 

A  second  approach  to  the  examination  of  applicable 
wave  theories  was  to  compare  measured  values  of  celerity  or 
wave  phase  speed  with  theoretical  values. 
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B.   EXAMINATION  OF  LOW  FREQUENCY  SPECTRAL  COMPONENTS 

Figures  (8)  and  (9)  give  the  results  of  the  spectral 
analysis  when  comparing  the  outer  and  center  profiles  and 
the  center  and  inner  profiles.   Figure  (10)  is  a  spectral 
comparison  between  horizontal  and  vertical  measured  ve- 
locity spectra.   Figure  (11)  is  a  superposition  of  the  three 
surface  profile  spectra.   Table  (3)  is  a  summary  of  low 
frequency  spectral  peaks  found  throughout  the  five  spectra 
generated  and  gives  a  wave  height  or  velocity  value  for 
each  peak  relative  to  a  maximum  spectral  peak  of  1.0.   These 
values  are  presented  in  this  manner  for  comparison  purposes 
only. 

The  following  observations  can  be  made  concerning  the 
results  shown  in  the  Figures  and  Table  (3) : 

1.  The  maximum  wave  energy  for  all  spectra  is  concen- 
trated at  a  primary  frequency  of  0.125  Hz  corresponding 

to  a  period  of  8.0  seconds. 

2.  A  secondary  peak  exists  at  0.098  Hz  (10.2  seconds). 

3.  An  approximate  subharmonic  of  the  primary  frequency 
exists  at  0.060  Hz  (16.7  seconds). 

4.  A  fourth  peak  of  interest  shows  up  at  0.011  Hz 
(92  seconds)  on  the  inner  profile  and  vertical  velocity 
spectra. 

5.  A  major  peak  on  the  horizontal  velocity  spectrum 
only  exists  at  0.016  Hz  (61.2  seconds). 
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SPECTRAw  -  OUTER  SURFACE  PROFILE  VS.  CENTER  SURFACE 
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6.   From  a  qualitative  examination  of  the  high  fre- 
quency ends  of  the  three   profile  spectra  (Figure  (11)),  it 
appears  that  turbulence  increases  to  a  peak  at  the  breakers 
and  the  waves  seem  to  reorganize  slightly  inside  this  point. 

Without  further  data  analysis,  discussion  of  these  ob- 
servations can  be  made  only  in  a  general  sense.   The 
observation  of  a  primary  frequency  at  0.125  Hz  agrees  well 
with  the  wave  periods  of  eight  to  ten  seconds  observed 
during  the  beach  experiment.   The  additional  peak  at  0.098 
Hz  indicates  that  incoming  waves  were  probably  composed  of  a 
combination  of  8.0  and  10.2  second  waves. 

Guzza  and  Davis  (1974)  state  that  in  a  theoretical 
analysis  of  edge  waves  excited  by  incoming  waves,  the  prom- 
inent edge  wave  mode  is  the  first  subharmonic  of  the  primary 
frequency.   An  energy  peak  was  observed  at  this  frequency. 
The  author  feels  that  the  peak  at  0.011  Hz,  which  shows  on 
the  inner  profile  and  vertical  velocity  spectra  only,  may 
also  be  caused  by  an  edge  wave  in  the  swash  zone.   Further 
experiments  with  profile  measurements  in  a  line  parallel  to 
the  beach  will  clarify  edge  wave  observations. 

The  comment  that  the  turbulence  inside  the  breaker  zone 
decreases  must  be  tempered  by  the  fact  that  the  inner  wave 
gage  was  not  the  same  type  as  the  other  two  gages.   It  is 
felt  that  the  Baylor  wave  gage  was  not  as  responsive  to  high 
frequencies  as  were  the  capacitance  wave  gages.   Future 
studies  with  matched  gages  across  the  entire  surf  zone  will 
clarify  this  question  of  wave  organization. 


37 


C.   EXAMINATION  OF  STATIONARITY 

As  the  tide  rises  and  the  mean  water  depth  increases, 
the  breaker  zone  should  move  shoreward,  and  the  mean  para- 
meters measured  (mean  depth  and  variance)  at  the  towers 
should  change.   This  series  of  events  would  result  in  a 
non-stationary  situation.   In  Section  (III. A)  it  was  stated 
that  if  the  data  were  ergotic,  the  measurements  taken  and 
analyzed  during  this  experiment  were  typical  of  all  waves 
in  the  breaker  zone.   This  rather  bold  assumption  can  be 
examined  by  checking  for  stationarity  of  the  data  over  the 
period  of  measurements. 

The  data  is  assumed  to  be  weakly  stationary  if  the  mean 
and  mean  square  values  calculated  are  constant  with  time  or 
if  the  auto-correlation  function  or  its  Fourier  transform, 
the  energy  density  spectrum,  is  constant  for  adjacent 
lengths  of  records.   Both  approaches  were  taken  to  verify 
stationarity.   Energy  density  spectra  were  computed  for  the 
outer  and  center  profiles  during  two  independent  30  minute 
time  segments  and  are  illustrated  in  Figure  (12) .   In  a 
general  sense  (examining  slopes  and  energy  peaks) ,  they 
compare  favorably. 

Mean  depths  and  variances  were  then  calculated  for  ten 
overlapping  samples  of  30  minutes  each  and  are  plotted  in 
Figure  (13) .   A  non-parametric  run  test  was  applied  to  the 
data  as  described  by  Bendat  and  Piersol  (1971)  to  check  for 
stationarity.   On  the  basis  of  this  test,  the  data  did  not 
qualify  as  stationary. 
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The  non-stationarity  can  be  explained  by  a  number  of 
observations.   First,  it  must  be  re-emphasized  that  the 
waves  measured  were  unusually  small  and  consistantly  of 
the  spilling  type.   The  bottom  hump  shown  in  the  beach 
profile  (Figure  (4))  at  -18m  appeared  to  trip  the  breakers 
in  the  vicinity  of  the  center  wave  gage  despite  an  apparent 
increase  in  the  still  water  depth.   The  height  of  waves 
outside  the  breaker  point,  which  can  be  related  to  the 
variance  of  the  outer  profile  measurements,  decreased  as 
the  depth  increased  with  time.   This  was  an  indication  of 
a  decrease  in  the  shoaling  effect  on  the  waves.   Set  down, 
the  difference  between  still  water  level  and  mean  depth  at 
the  breakers,  increased  as  the  breaker  wave  height  decreased 
at  the  center  wave  gage.   Future  analysis  should  concentrate 
on  this  set  down  effect  and  its  relation  to  set  down  theory. 

Thus,  the  rising  of  the  tide  was  the  main  factor  contri- 
buting to  the  non-stationarity  of  the  data.   Despite  the 
fact  that  the  data  was  non-stationary,  the  analysis  was  still 
valid  when  the  fact  that  the  data  applies  to  a  very  special 
set  of  circumstances  is  considered. 
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D.   LINEAR  THEORY  APPLICATIONS 

Airy  or  linear  wave  theory  assumes  that  the  surface 
profile  is  sinusoidal,  i.e., 

n(t)  =  A  cos  (at)  (1) 

Linear  theory  horizontal  and  vertical  water  particle 
velocities  at  the  depth  z  (z  positive  upward)   below  the 
mean  water  level  are  given  by 

,/*.\    *~  cosh  (k(h+z))      .„..  #_, 

U(t)  =  Aa   sinh  (kh)    cos  (CTt)  <2> 

and 

w(t)    =  Aa   Sinh  ik^fn    sin    (at)  (3) 

sinh    (kh) 

Horizontal   and  vertical   velocity   spectra  can  be  computed 
from  the  wave   spectrum  using   the   equations 

q      (f)   -     |~cosh    (Mh+z))1       c    (f)  (4) 

Suc(f)    -     L      sinh    (kh) J        SnUJ  (4) 

and 

q  (f)    -      fsinh  (k(h+z))l   „  (f)  f5, 

Swc(f)  -   L   sinh  (kh) J   S^t]  (5) 

The  wave  number,  k,  was  calculated  using  an  iteration  of 
the  linear  theory  characteristic  equation 

a2  =  gk  tanh  (kh)  (6) 
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Calculated  horizontal  (S   )  and  vertical  (Swc)  velocity 
spectra  were  computed  using  the  center  surface  profile 
spectrum  and  compared  with  measured  spectra  (S,^  and  Swm) 
as  shown  in  Figures  (14)  and  (15).   The  calculated  and 
measured  horizontal  velocity  spectra,  Suc  and  S^^,     (Figure 
(14))  compare  well,  indicating  that  linear  wave  theory  was 
applicable  for  this  experiment.   On  the  other  hand,  the 
calculated  and  measured  vertical  velocity  spectra  (Swc  and 
swm^ '  shown  in  Figure  (15)  did  not  agree  well  at  all.   The 
lack  of  agreement  between  these  vertical  velocity  spectral 
magnitudes  could  be  partially  attributed  to  a  low  signal  to 
noise  ratio  in  the  measurements  (note  the  vertical  velocity 
time  sample  relative  to  horizontal  velocity  time  sample  in 
Figure  (5)  which  greatly  reduced  the  accuracy  of  the  measured 
velocity  spectrum. 

Phase  measurements  are  used  to  substantiate  the  findings 
that  linear  theory  was  applicable.   Equations  (1)  and  (2)  show 
that  the  center  surface  profile  and  horizontal  velocity 
measurements  should  have  been  in  phase  (both  have  cos  (at) 
terms) .   An  approximately  zero  degree  phase  was  observed  well 
beyond  the  0.5  coherence  point  (Figure  (14)).   The  cos  (at) 
and  sin  (at)  terms  of  equations  (1)  and  (3)  indicate  that  the 
surface  profile  and  vertical  velocity  measurements  should  be 
90  degrees  out  of  phase.   Figure  (15)  illustrates  that  this 
phase  difference  was  observed  in  a  range  up  to  the  0.5 
coherence. 
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Finally,  a  high  coherence  would  indicate  a  high  degree 
of  linearity  in  the  transfer  functions  between  different 
spectra.   Coherence  above  0.5  was  observed  for  an  approxi- 
mate frequency  range  of  0.0  to  0.65  Hz  (Figures  (14)  and 
(15) )  .   This  high  value  again  substantiated  the  excellent 
results  obtained  when  a  linear  wave  theory  transfer 
function  was  applied  to  the  center  surface  profile  spectrum 
when  calculating  the  horizontal  velocity.   Whereas  coherence 
is  high  for  the  same  range,  linear  wave  theory  grossly 
overpredicts  the  vertical  velocity  spectrum.   The  author 
is  led  to  believe  that  a  linear  relationship  between  the 
surface  profile  and  vertical  velocity  exists  but  is  not 
predicted  well  by  linear  wave  theory.   A  further  indication 
of  the  linearity  between  measurements  was  noted  due  to  the 
high  coherence  between  horizontal  and  vertical  velocity 
measurements  as  shown  in  Figure  (10) . 

E.   CELERITY  COMPARISONS 

A  summary  of  the  various  celerity  calculations  made  in 
this  section  is  included  in  Table  (4) .   The  measured 
celerity,  C, ,  is  computed  considering  the  shift  in  phase 
existing  between  two  adjacent  wave  gages  as  plotted  in 
Figure  (16) . 
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At  x-l  =  0,  let 

nl  =  Al  cos  (~at)  (7) 

Then  at  x„  =  x 

H2  =  A2  cos  (kx-crt)  (8) 

Therefore,  the  phase  shift,  <j> ,  between  stations  1  and  2 

becomes 

(J>   =  k(f)    x  (9) 

By  definition 

C(f)    =    <L£J>  (10) 

k(f ) 

Or  substituting  and  solving  for  C  (f) 

C  (f )    =    °  ^f )  x  (11) 

The  dotted  lines  in  Figure  (16)  plot  a  constant  change  in 
phase  with  frequency  over  the  coherent  ranges  of  Figures 

(8)  and  (9) .   This  surprising  result  indicates  that  wave 
celerity  in  the  surf  zone  did  not  depend  on  frequency,  that 
is,  the  waves  were  non-dispersive. 

The  celerity  was  then  computed  for  points  between 
stations  using  the  linear  theory  approximation  for  shallow 
water  waves 

C2  =  (gh)^  =  constant  (12) 

These  values  agree  well  with  the  measured  celerity,  C. 

(+  4  percent).   Again,  as  stated  in  the  previous  section,  it 
appears  linear  theory  is  applicable  in  the  surf  zone. 
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Using  the  significant  wave  height,  H  ,  calculated  from 

s 

the  variance,  celerity  was  computed  using  the  approximation 
for  both  cnoidal  theory  and  solitary  wave  theory 

C3  =  (g(h+Hs))^  (13) 

These  values  overpredict  the  measured  celerity  by  20  percent, 

As  a  final  step,  celerity  was  calculated  as  a  function 
of  frequency  using  equation  (10)  for  linear  wave  theory. 
The  wave  number,  k,  was  calculated  from  an  iteration  of  the 
linear  theory  characteristic  equation,  equation  (6) .   The 
resulting  celerities,  C. (f ) ,  are  graphed  in  Figure  (17). 
The  decrease  in  celerity  with  increasing  frequency,  indica- 
tive of  a  dispersive  wave  system  and  expected  with  linear 
theory,  would  result  in  a  corresponding  increase  in  phase 
shift  with  frequency.   This  relationship  was  not  observed; 
rather,  the  phase  shift  with  frequency  was  constant  as 
illustrated  in  Figure  (16) .   The  value  for  C4  (f )  at  the 
primary  frequency,  0.125  Hz,  agreed  well  with  the  measured 
celerity,  Ci  (+  2  percent) . 

In  summary,  celerity  calculations  created  a  dilemma 
which  future  studies  in  the  surf  zone  will  have  to  resolve. 
The  constant  phase  shift  and  measured  celerity  indicated 
that  all  frequency  components  are  non-dispersive  and  ride 
on  the  back  of  a  non-linear  wave  form  moving  at  the  speed 
C-,.   However,  linear  wave  theory  predicted  the  horizontal 
velocity  and  both  C2  and  C4  (0.125Hz)  closely  approximate 
C^ .   If  linear  theory  is  applicable,  the  wave  system  is 
dispersive  or  each  component  moves  at  a  celerity  which 
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depends  on  its  frequency.   It  appears  that  while  linear 
theory  can  predict  the  relationships  between  waves  it  cannot 
be  used  to  define  the  breaker  zone  wave  form  which  is  non- 
linear and  non-dispersive. 

F.   COMMENTS  ON  INSTRUMENTS  AND  BEACH  EXPERIMENT 

As  a  corollary  goal  for  this  experiment,  the  instrumen- 
tation used  and  its  placement  in  the  surf  zone  was  to  be 
tested.   The  capacitance  wave  gages  used  were  found  to  be 
reliable,  inexpensive,  hardy,  and  easily  monitored  in  the 
surf  zone.   They  were  very  accurate  and  appeared  to  give  an 
excellent  representation  of  wave  profiles  as  they  actually 
occurred.   The  work  required  to  construct,  calibrate,  and 
install  the  gages  was  minimal.   By  mounting  all  electronics 
in  a  sealed  container  and  running  signal  cables  overhead,  the 
chance  of  shorting  the  system  and  introducing  noise  was 
greatly  reduced. 

The  use  of  electromagnetic  flow  meters  in  the  surf  zone 
gave  excellent  results.   The  lack  of  moving  parts  which  clog 
with  sand  and  the  ability  to  record  and  monitor  signals  on 
the  beach  virtually  eliminated  the  problem  encountered  by 
previous  experimenters  of  incomplete  or  questionable  records. 
Calibration  of  the  flow  meters  was  somewhat  difficult  and 
time  consuming  and  required  special  equipment.   A  new  cali- 
bration scheme  should  be  devised  to  increase  the  calibration 
frequency  range  and  accuracy.   The  results  obtained  were 
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sufficient  for  the  experiment,  however.   Calibration  of 
the  flow  meters  is  required  as  their  response,  although 
linear,  deviated  from  the  manufacturers  specifications, 
showing  both  a  slight  DC  bias  and  a  different  calibration 
slope. 

The  Baylor  resistance  wave  gage  also  proved  reliable, 
accurate,  and  very  easy  to  calibrate.   On  the  beach  it  was 
unwieldy  and  required  a  special  tower  to  maintain  tension  on 
the  cables.   It  was  felt  that  the  Baylor  gage  output  was  not 
as  accurate  as  that  obtained  with  the  capacitance  gages. 

The  selection  of  instrument  locations  on  the  beach  was 
conducted  four  days  prior  to  the  actual  measurements  and 
was  done  by  considering  the  beach  profile  and  the  location 
of  breakers  on  that  day.   Wave  gages  were  placed  to  span 
the  breakers  with  the  velocity  measurements  to  be  taken  in 
the  breaking  waves.   Fortunately,  wave  characteristics 
remained  constant  until  the  date  of  the  experiment  and 
placement  was  correct  for  the  intended  measurements . 

The  key  word  for  the  experiment  setup  was  organization. 
The  equipment  and  manpower  requirements  for  an  examination 
of  surf  zone  kinematics  using  temporary  stations  can 
quickly  become  overwhelming.   It  is  felt  by  the  author  that 
five  to  six  measurements  are  maximum  for  a  one-man  experiment 
in  the  surf  zone. 
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V.   CONCLUSIONS 

With  the  conditions  of  a  mild  surf  and  spilling  breakers 
linear  wave  theory  was  shown  to  predict  wave  kinematics  well. 
Spectral  analysis  showed  a  good  agreement  between  the 
measured  horizontal  water  particle  velocity  spectrum  and 
the  horizontal  velocity  spectrum  calculated  using  a  linear 
theory  conversion  factor  on  the  surface  profile  spectrum 
over  the  highly  coherent  band  of  prominent  wave  energy.   A 
similar  high  coherence  and  the  agreement  with  predicted 
phase  relationships  over  the  same  band  found  in  a  cross- 
spectral  analysis  of  the  surface  profile  and  vertical  water 
particle  velocity  led  the  author  to  believe  that  a  linear 
relationship,  although  not  predictable  by  linear  wave 
theory,  also  existed  which  relates  the  vertical  velocity 
spectrum   to  the  surface  profile  spectrum. 

The  shallow  water  approximation  for  celerity  used  in 
linear  theory  was  an  excellent  prediction  for  wave  phase 
speed  in  the  experiment  environment.   A  constant  phase 
shift  with  frequency  indicated  the  breaking  waves  were 
non-dispersive  or  that  celerity  was  not  frequency  dependent. 
Apparently  all  energy  rode  on  the  back  of  a  non-dispersive, 
and  therefore  non-linear,  waveform  which  was  at  the  primary 
frequency.  An  inconsistancy  thus  arose  in  that  linear  theory, 
a  dispersive  theory,  could  be  used  to  predict  wave  kine- 
matics; yet,  celerity  calculations  all  pointed  to  a 
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non-dispersive  wave  form.   Future  studies  in  the  breaker 
zone  are  needed  to  resolve  the  question  of  linearity. 

The  use  of  capacitance  wave  gages  and  electromagnetic 
flow  meters  is  a  convenient  method  of  gathering  accurate, 
continuous  data  in  the  surf  zone. 
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